Instabilities and self-oscillations in atomic four-wave mixing 
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The development of integrated, waveguide-based atom optical devices requires a thorough under- 
standing of nonlinear matter-wave mixing processes in confined geometries. This paper analyzes the 
stability of counterpropagating two-component Bose-Einstein condensates in such a geometry. The 
steady state field equations of this system are solved analytically, predicting a multivalued relation 
between the input and output field intensities. The spatio-temporal linear stability of these solutions 
is investigated numerically, leading to the prediction of a self-oscillation threshold that can be ex- 
pressed in terms of a matter-wave analog of the Fresnel number in optics. 
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I. INTRODUCTION 



The recent development of narrow atomic waveguides 
micro-fabricated on glass chips [Q-Q has raised the ex- 
citing possibility of the design and manufacture of inte- 
grated atom-interferometry-based sensing devices. With 
the inclusion of an 'atom laser' (5[|| as a high-brightness 
source of coherent atomic matter-waves, it is possible to 
imagine 'practical' devices which could compete with or 
out-perform conventional optical interferometric sensors. 
The use of high-density atomic fields comes at a price, 
however, as atomic matter waves are subject to nonlin- 
ear wave mixing due to atom- atom interactions. It is of 
crucial importance, therefore, to understand the effects of 
nonlinear wave mixing on wave-guide based atom-optics 
devices so that they may eventually be controlled or even 
exploited. In this spirit, the present paper is a first at- 
tempt at an analysis of wave-mixing instabilities in quasi- 
one-dimensional ultracold atomic samples. 

The observation of atomic four-wave mixing (^-^] 
and solitons |Tfj|-|20|] in Bose-Einstein condensates of di- 
lute atomic vapors pl}-|23(] clearly demonstrate both the 
significance of nonlinear effects in quantum degenerate 
atomic fields, as well the benefits of exploiting the mathe- 
matical analogy between the nonlinear equations describ- 
ing self-interacting Schrodinger fields and those describ- 
ing the propagation of light in nonlinear media. Non- 
linear wave-mixing instabilities have been studied exten- 
sively in nonlinear optics, and many of the techniques and 
results developed can readily be adapted to the problem 
at hand. 

Focusing on effective one-dimensional geometries, the 
question of stable and unstable steady-state configura- 
tions has long been a topic of optical research. Win- 
ful and Marburger |24j first proposed that bistability 
could occur in collinear degenerate four-wave mixing and 
shortly thereafter Silberberg and Bar- Joseph 1 25 1 showed 
that even for the rather simple case of equally polarized 
counterpropagating laser beams instabilities and even 
chaos may occur in the dynamical behaviour. Multi- 
branched steady-state solutions were first derived by Ka- 
plan and Law pa, however the stability of the steady- 



state field configurations was not determined. Consid- 
erable work on the spatial or temporal stability of such 
systems was subsequently carried out by many others 
p7|-p2[. In particular, optical instabilities and polariza- 
tion bistability were experimentally observed in sodium 
vapor [p||3T 
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FIG. 1. Four matter waves incident into a region of non- 
linear interaction. The two forward moving and the two back- 
ward propagating modes, with opposite wave vectors, are dis- 
tinguished by their internal state. 

Bistability and nonlinear instability are typically re- 
lated to four-wave-mixing phenomena in systems which 
exhibit a cubic nonlinearity. In ultracold atomic sys- 
tems it is readily shown that in the s-wave scattering 
approximation the form of the self-interaction is that of 
a cubic nonlinearity. A recent paper by Law and cowork- 
ers [ |35| analyzed four- wave-mixing processes between the 
hyperfine ground-state components of components of a 
23 Na spinor condensate confined in an optical dipole trap. 
Goldstein and Meystre |3(| presented a full quantum- 
mechanical theory of four- wave mixing in a system where 
two mp = momentum side modes were counterpropa- 
gating while the mp = ±1 states were at rest. 

In the present work we investigate a collinear four- 
wave mixing geometry as sketched in Fig. [j], where each 
of the counterpropagating matter waves can be in one 
of two different atomic states, for example two hyperfine 
state levels of 87 Rb. This situation is closely analogous 
to the optical case, the two internal atomic states taking 
the place of the polarizations of the field. As such, this 
system is formally equivalent to the case of counterprop- 
agating light fields in an Kerr medium. 

In contrast to the exact quantum treatment of 
Refs. [j3oj|36|, our analysis is based on a mean- field ap- 
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proach, the matter-wave equivalent of treating the elec- 
tromagnetic field classically. We investigate both the 
steady-state and dynamical behavior of the system by 
a combination of analytical and numerical methods. The 
output fields are found to generally exhibit a multi- valued 
dependence on the inputs, characteristic of bistable and 
multi-stable systems. The stability analysis for this par- 
ticular configuration, however, shows that only the upper 
branch of the steady-state curve is stable against small 
perturbations. More interesting, however, is the occur- 
rence of a threshold behavior in the output fields indi- 
cating the onset of self-oscillations in the system. The 
feedback mechanism leading to this effect is the grating 
established in the medium by the interference between 
the various fields. 

We note that matter- wave bistability was recently pre- 
dicted in a simple model of a driven nonlinear Gross- 
Pitaevskii equation J3?J which neglected however the ef- 
fects of collisions between the strong driving field and 
the condensate. In contrast, the present system includes 
both the effect of two-body collisions, and fully accounts 
for propagation effects. 

This paper is organized as follows: Section || intro- 
duces our model and derives the nonlinear partial differ- 
ential equations describing the propagation of the inter- 
acting atomic beams. Section 111 solves these equations 
analytically in steady state and shows the appearance of 
multistable solutions and threshold behavior. The stabil- 
ity of the steady state solutions is investigated in Section 
IV using numerical methods, while section M discusses 
the onset of self-oscillations. Finally, Section VI is a sum- 
mary and conclusion. 



II. MODEL 

We consider an ultracold two-component Schrodinger 
field *(r,t) = (*i(r, t), # 2 (r, t)) T , the indices 1 and 2 
labeling the internal state of the atoms of mass to, e.g. 
two hyperfine ground states. These fields consist of two 
counterpropagating plane waves which interact propagat- 
ing along the z-axis in an interaction region < z < L, 
the nonlinear interaction outside this region being turned 
off e.g. by tuning a magnetic field to a Feshbach reso- 
nance. 

Our starting point for the description of this system 
is the many-body Hamiltonian describing the evolution 
of a two-component condensate, the effects of collisions 
being described in the s-wave scattering approximation, 



H=J2 M 3 ^jM)(f- + ^)^(r,i) 



tiff tiff 



2g x ¥ 1 ¥ 2 y 1 i> 2 



(1) 



where the scattering strengths gi are related to their re- 
spective s-wave scattering lengths at by 



Antia; 



fji 



(2) 



For T — > 0, the condensate is well described by 
a two-component Hartree condensate wave function 
4>{v,t) — (0i (r, t), 02( r i t)) T : governed by the coupled 
Gross-Pitaevskii nonlinear Schrodinger equations 



ih$! = ( -—V 2 + Vi ) $i 

2m 



-HN [g s \$i\ 2 + g x \<h\ 2 ] $i 



(3) 



and similarly for $2 with 1 <-» 2, where the total atomic 
density is normalized to 1 by factoring out the number of 
atoms N. Quantum fluctuations about this solution can 
be analyzed by introducing the mean-field approximation 



*i(r,t)~$i(r, t)+5^(r,t), 



(4) 



where the bosonic operator Si/j^r, t) describes small fluc- 
tuations about the mean field $j(r,t) = (^(r, t)). This 
analysis will be the subject of future work. 

We assume that the atomic fields are tightly confined 
in the transverse dimension but free to move in the third 
one such that the motional degrees of freedom in the 
x — y plane are frozen, a situation that could be realized 
in atomic waveguides. In that case, we may factorize the 
ground state Hartree wave function into a parallel and a 
transverse part as 



$ j (r,t) = <f> { l ) (x,y)<l> j (z,t)e- 



(5) 



where (j)^{x,y) is taken to be the normalized ground 
state of the transverse potential Vj{x,y) with energy 
fwjj. The problem is then reduced to an effective one- 
dimensional geometry, with the longitudinal condensate 
wave function satisfying the one-dimensional nonlinear 
Schrodinger equation 



ih(pi(z, t) 



h z d 2 



4>i(z,t) 



2m dz 2 
+HN [ mgi \M^t)\ 2 + 

and similarly for 4>2(z) with 1 <-» 2. Here 
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dx dy 
dx dy 



.7 = 1,2 
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(7) 



We consider the situation where two counterpropagating 
beams of matter waves are moving along the axis of the 
waveguide, with 



<f>j(z, t) = [<f>jf(z, t)e lkz + <f> jb (z, t)e 



—ikz 



(8) 
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where j = 1,2 and hk 2 /2m = u>. We assume that the 
spatial envelopes of these beams vary slowly over a de 
Broglie wavelength, the atom optics version of the slowly 
varying envelope approximation, 



(9) 
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1/J 



(11) 



as well as two additional equations with 1 <-» 2. In 
these equations we have scaled the nonlinear coupling 
constants to the overlap integral g and the total particle 
number N as 



9j 
9x 



VjNgj, j 
VxNg x . 



1,2 



(12) 



The factors of 2 appearing in Eqs. ( |Io| ) and (|TT|) result 
from the nonlinear nonreciprocity familiar in nonlinear 
optics. 



III. STEADY STATE 

To find the steady state of the system, we proceed by 
first observing that the total atomic density 



Q : 



9 i/l 



'16 



+ |02/| 2 + |026| 5 



(13) 



is a constant of motion. For simplicity we take g\ — 92 — 
g x = g (this can be achieved by tuning the scattering 
lengths via e.g., Feshbach resonance and/or by adjusting 
the transverse potential Vy); in the optical analog this 
corresponds to a purely electrostrictive Kerr medium. 
Dropping the time derivatives from Eqs. ( p~0| ) and ( pd| ) 
and introducing the scaled velocity 



gives then 



dz 



>ib 



dz 



hk 
mg 



il / + R(j)ib 
!>16 + i?*0i/, 



(14) 



(15) 



and 



d(j)2f 

iv— — = g<p 2 f + Rq>2b 



dz 

#26 

dz 



g4>2b + R*4>2f, 



(16) 



where the introduction of the new variable 

R(z) = <tn f <l>Z b + <hf<t>$b (1 7 ) 

allows us to decouple the evolution of the two internal 
components of the field two condensate species 1 and 2. 
This can be seen from the observation that 



(10) ^ = ^il (b *,, f ^u + dh L(h * 



dz dz 



lb I" 91 f 



dz dz 



•>2b^92f 



d0* 



21, 



dz 



-3gR 

v 



which yields 



R(z) = i?oe~" 3<?z . 



(18) 



(19) 



We can therefore determine the steady state for the two 
internal states of the condensate separately. 
Substituting the new rotating field amplitudes 



ipif = <j) 1} e- ez 
tpib = 0i&e~" ez , 



(20) 



into Eqs. (|15|) yields the second-order differential equa- 
tion 



d 2 ip lf igdipif \R \ 2 

=:: u —l/j lf , 



dz 



v dz 



(21) 



which shows that the propagation of ip± / is characterized 
by the spatial frequencies 



Q , 1 



It can be shown that 



M 2 = g 2 -4|i? | 2 > 



(22) 



(23) 



so that ipif(z): 



, M \ (M 
Aif sm [ — z J + cos ! — : 



2v 



(24) 



the constants A\ / and B\ / being determined by the 
boundary conditions. We observe that A\ f depends on 
both ip±b (0) an d Rq, and hence on the boundary condi- 
tions for all four fields, 0i/(O), <fiib(L), 02/(0) and 2 &(£)- 
Therefore the equations of motion of all four fields need to 
be solved and used to calculate the respective coefficients 
for any one field. The explicit forms of the coefficients 
A^i and B U j, where fi — 1,2 and % = f,b are given in 
appendix [A]. 
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Further analytical progress can be achieved by decom- 
posing the field ipif into a real amplitude and phase fol- 
lowing reference p6fl, 



^1/ = V^i7 ex p(^i/)> 



(25) 



and concentrating on the field amplitudes only. One finds 
readily 



Pif( z ) 



2 -( Pl/ (0)-/3) z sin 
M 



+ (pi/(0) -0) cos 



M 



-z 



r 



(26) 



where 



Bi f ) 



\Bi f \ 2 ), (27) 

and the sign in front of the square root in Eq. (^) is 
determined by the sign of AifB^ + A*jBif. Similar 
relations hold for the other field components. 

One advantage of concentrating on the amplitudes p^i 
only is that it is sufficient to know one of them to deter- 
mine the others. This follows from the fact that 



dpif _ dpib _ dp 2 f 



dz 



dz 



dz 



dp2b 
dz 



(28) 



which allows us to introduce the three conserved quanti- 
ties 



PS = PiS( z ) + Pvi z ) 
Pb = Pib(z) + P2b(z) 
p x = p lf (z) + p 2b (z), 



(29) 



so that 



Pib(z) 

P2S(Z) 
P2b(z) 



Pb 

PS 

Px 



p x + pif(z) 
- Pif{z) 
■ Pif{z)- 



(30) 



The existence of these conservation laws was previously 
poi nted out in the context on nonlinear optics in Refs. 
@ and @. 

For concreteness, we now consider the specific exam- 
ple where the intensities of the forward and backward 
propagating fields are equal, 



PS = Pb, 



and 



Pib(L) 

P2b(L) 



Pb 
0. 



g/2 



(31) 



(32) 



Under these conditions, one finds readily that 



and 



\Ro\ 2 = PbPif(L) = -gpif(L) 



\a\=p = p lf (L)/2. 



(33) 



(34) 



With Eqs. (g3|) and (|3J), Eq. (|26j) reduces to the remark- 
ably simple form 



where 



Pi/(0) 
PisW 

{ 2mg\ 



cos 2 (kL/2) 



(35) 



« = M/v = (j^j g^l/4-p lf (L)/2g. (36) 

1/k defines a characteristic length. Equation ( |35| ) pre- 
dicts a multivalued relationship between the input and 
output intensities p\ f . The longer k, the higher order of 
"multistability" predicted by Eq. (p5|). Fig. 2 illustrates 
the input-output relationship of pif. 

Alternatively, and recalling Eqs. ( |l2] ) and ( |l4| ) we ob- 
serve that the argument kL/2 of the cosine function in 
Eq. ( |35| ) is also proportional to the total density of atoms. 
This leads us to expect some interesting behaviour when 
varying the number of particles involved, a scheme whose 
optical analog has been the object of considerable work 
p6| , p8| -|3C|] . However, we don't further investigate these 
questions in the present paper, which is limited to the 
case of fixed total atomic density. 
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FIG. 2. Normalized output intensity p±f(L) versus normal- 
ized input intensity pi/(0) from Eq. (Ba). 1: L/L c = 0, 
2: L/L c = 2tt, 3: L/L c = 4tt and 4: L/L c = 8tt, with L c 
as defined in {J39|) . 



Returning from the scaled version (|2[ ) of the coupling 
constant g to its original definition (g) via Eqs. (0), we 
observe that for a transversely homogeneous sample, the 
factor gg in Eq. (|3^) becomes ggv, where gy is the vol- 
umetric density of the atomic system (as opposed to its 
linear density g). Eq. (5^) becomes then 



4 



= k 



2m 
Wk* 
E m f 



hg gv \ l/4-2pif(L)/2g 



yJl/A- p lf {L)/2 e . (37) 



where we have introduced the kinetic energy Ek e — 
h 2 k 2 /2m and the mean-field energy E m f = hggv- 
Introducing the healing length 



2mgg v 



(38) 



o 




and the de Broglie wavelength X^b — 2n/k shows that the 
multistable properties of the system are fully determined 
by the characteristic length 



E, 



kE, 



4tt j 



mf 



\lb 



(39) 



IV. STABILITY ANALYSIS 



In this section, we analyze the stability of the steady- 
state solution (|3^) against small classical perturbations. 
For simplicity, we assume that the fields at the bound- 
aries of the interaction region are real, 



hf{0) = yfpTj 

Mo) = Vp7 

kb{L) = 

f>2b(L) = yfpb- 



Plf 



(40) 



FIG. 3. log (MA r)) Pi fss (L) = 0.2/9/. S lf (L,r) is 
scaled to the intensity p/ of the forward moving field and the 
time t is in units of the round-trip time through the interac- 
tion region. Curve (a) is an example of unstable behaviour at 
L/L c = 4-7T (nL/2 ~ 0.897r), while curve (b) shows stability 
at L/Lc = 1-67T (kL/2 x 0.36tt). 

Both approaches indicate that the only stable branch 
of the solution ( |35| ) is the "uppermost branch" , i. e. the 
branch corresponding to the highest pif(L) for a given 
pi / (0) which corresponds to the condition 



kL 

< — < tt/2. 



(42) 



Figure ^ shows a typical result for the temporal evo- 
lution of 8if at z — L in both an unstable and a stable 
branch. After some short-time transients, the unstable 
dynamics becomes completely dominated by the largest 
positive eigenvalue and the growth of the small pertur- 
bation becomes exponential. 



V. SELF-OSCILLATIONS 



The first approach proceeds by expressing the conden- 
sate wave function <pif(z,t) and its complex conjugate 
as 

cj> lf (z, t) = 4> lfss (z, t) + S lf {z)e~ lXt 
ct>i f (z,t)=<t>* lfss (z,t) + e lf (z)e- iXt (41) 

and linearizing the equations of motion (|l0|) and ( |IT| ) , and 
similarly for the other field components. The spectrum of 
eigenvalues A can be found by discretizing the resulting 
system of equations with N points along the z-axis which 
yields a sparse 8iV x 8N eigenvalue problem that can 
be solved numerically using standard techniques. The 
steady-state solution is then unstable only if eigenvalues 
with positive imaginary parts appear in the spectrum. 

In addition to this temporal stability, we have also car- 
ried out a full spatio-temporal analysis by solving directly 
the linearized form of Eqs. ( |l0| ) and (|ll]) for small per- 
turbations about the steady state of the fields. 



An important consequence of the stability analysis is 
the prediction of a self-oscillation threshold in the system, 
as can be seen by considering the case pi/(0) = 0, that 
is, no input field in mode "1/". In that case, Eq. (|35| ) 
reduces to 

Pi/(L)cos 2 (kL/2) = 0. (43) 

Hence in the stable region given by ( f42| ) , p\ / (L) must be 
0; while in the unstable region, pif(Lyma.y take finite val- 
ues which results from the amplification of fluctuations. 
The onset of instability is given by 

kL = kL ^l/4- Plf (L)/2g = tt. (44) 

Since the argument of the square root is always less than 
1/4, we have that 
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Hence the threshold condition (H) cannot be met unless 



E, 



mf 



Ei.p 



> 



2tt 
kL' 



i\ 1 
8tt 3 



LX. 



lib 



which can be recast as 



L 



> 2?r. 



(46) 



(47) 



(48) 



If this threshold condition is satisfied, the field <f>\f 
undergoes a second-order-like phase transition to self- 
oscillations, as illustrated in Fig. [|. 

The physical origin of the self-oscillations is the dis- 
tributed feedback resulting from the cross-phase modu- 
lation with the other fields present. In order to lead to 
gain at the wavelength of 4>if, this grating must itself 
be modulated with period 2n/k. But the mean-field en- 
ergy of the condensate fights against the creation of such 
spatial modulation. The healing length, whose associated 
momentum i^ 1 yields a kinetic energy contribution equal 
to the mean-field energy, is the smallest length scale over 
which the required changes can occur. Hence, it is not 
surprising that the threshold condition should be related 
to the healing length, the length of the sample and the 
atomic de Broglie wavelength. 



situation is different in that we now consider the longitu- 
dinal stability of the system. Still, the analogy is rather 
telling. For samples of length short or comparable to the 
healing length, the condensate is so stiff as to prevent 
the generation of higher longitudinal modes, hence the 
instability requires a condensate much larger than t^. As 
such, i^/XdbL can be thought of as the longitudinal Fres- 
nel number of the condensate, with the healing length 
playing a role similar to that of the system aperture a in 
optics. 

The onset of self-oscillations is illustrated in Fig. [|, 
which shows the evolution of pif(L,t) for /9i/(0,t) = 
0.0 lg/ 2, this small value simulating some small fluctua- 
tion about pi/(0,i) = 0. This simulation assume that all 
field amplitudes inside the interaction region (0 < z < L) 
are initially equal to zero. 
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FIG. 5. pi f (L, t) for pi/(0) = O.Olg. The time is in units 
of the round-trip time through the interaction medium. Curve 
(a) is in the regime with only one high output stable branch 
at L/L c = 4-7T; curve (b) is at L/L c = 1.6-7T where for every 
(low) output intensity there exists a stable steady-state field 
configuration. 



VI. SUMMARY AND CONCLUSION 
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FIG. 4. Self-oscillation threshold of the intensity pi / (L) as 
a function of the interaction length L. Above threshold, this 
field spontaneously builds up from fluctuations even in the 
absence of input, pi/(0) = 0. 

It is interesting to note that the quantity t\/\dbL is 
reminiscent of the Fresnel number T — a 2 /XL in optics, 
where a is the aperture of the system, A the wavelength 
and L the distance of propagation. The Fresnel num- 
ber is a measure of the number of transverse modes that 
can be excited in an optical system. For very large Fres- 
nel numbers, the wave can well be approximated as a 
plane wave, while diffraction effects and multiple trans- 
verse modes become important for small T . The present 



In summary, we have studied a system of counterprop- 
agating two-component Bose-Einstein condensates in a 
waveguide configuration. This system exhibits interest- 
ing nonlinear behavior such as four-wave mixing and self- 
oscillations. We have presented an analytical solution of 
the steady state field equations and found a multivalued 
relationship between the input and the output field in- 
tensities. Through a linear perturbation analysis as well 
as a full numerical study, we have found that one and 
only one branch of solutions is stable for a given input 
configuration. The onset of the instability depends on a 
parameter in close analogy with the Fresnel number in 
optics. 

We remark that these results were obtained under the 
assumption that all the effective scattering lengths are 
the same (gi = g% = g x ). In future work, we plan to 
investigate the case of unequal scattering lengths, which 







will make our system formally equivalent to the case of 
counterpropagating light fields in a non-electrostrictive 
Kerr medium. The analogy with the optical system also 
suggests that it would be desirable to allow for different 
intensities of the two counterpropagating matter-wave 
fields. It has been shown that bi- or multistable solu- 
tions exist in these nonlinear optical systems. It will be 
interesting to see if such a matter wave system can also 
exhibit multi-stability. 

Finally, it will be worthwhile to extend the analysis 
past the Hartree mean-field theory in order to study the 
onset of nonclassical effects, such as the squeezing asso- 
ciated with the anomalous density of the system. 

The experimental realization of this system is certainly 
challenging, but we see no fundamental problems using 
present day technology. Indeed, the question of stabil- 
ity in wave-mixing processes will most likely arise natu- 
rally in the next generation of atom interferometry ex- 
periments. 



B 2 b = 



Mb = 



i2B*fBi b B 2 f 



,-3i e L/2v M( t>2b{L) 

sin(A/L/2u) 



9(\B2f\ 2 ) 
4>ib(L)e~ 3ieL / 2v — i?ibCOS (Jf-fc) 



A 2b = -B 2b /t&n(ML/2v) 

M f = ^ [gB lf - 2 (B lf B* b + B 2f B* 2b ) B lb ] 

A 2f = ^ [gB 2f - 2 (B lf B* lb + B 2f B* 2b ) B 2b ] (A2) 



where 



g(x) = i(g- 2x) - M/ tan(ML/2v). 



(A3) 



These equations simplify considerably in the specific case 
02b (L) = that we have analyzed in detail. Similar, but 
simpler equations are easily derived if one of the input 
fields is equal to zero or for sin(ML/2v) ^ 0. 
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APPENDIX A: 

Solving Eq. ( |2l| ) and its counterpart for the backward 
moving field and inserting the results into Eqs. (20) gives 



4>u{z) 



■ 3g 



M \ (M 
A 1/S m[ cob (-* 



■ 3o 



. M \ (M 
A u sin — z + B lb cos — z 
2v ) \2v 



(Al) 



and similarly for 1 <-> 2. The boundary conditions are 
^1/(0), 0i6 (L), 02/(0) and <f> 2b {L). Since the four fields 
are coupled by i?o = BifB* b + B 2 jB^ h it is necessary to 
consider all fields to determine the coefficients A^i and 



B^. We start from Eqs. (Al), we derive expressions for 
0j/(O), d<pjf(0)/dz, 4>jb{L) and d<pj b (0)/dz, j = 1, 2 in 
terms of these coefficients and use the relations ([H]) to 
close this set of equations. 

For sin(ML/2u), 0i/(O) and 02/ (0) ^ we find 



Bif = 01/(0) 

B2f = 02/(0) 

B\ b = — - 



V 9{\B 2 fr) 



shx(ML/2v) 



g(\B lf \ 2 ) 



4 B 



i/l 



B 



2/1 



g{\B2 f ?) 
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